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Abstract The main goal of the paper is to establish time semidiscrete and space-time fully discrete maximal 
parabolic regularity for the time discontinuous Galerkin solution of linear parabolic equations. Such estimates 
have many applications. They are essential, for example, in establishing optimal a priori error estimates in non- 
Hilbertian norms without unnatural coupling of spatial mesh sizes with time steps. 

Keywords maximal parabolic regularity • finite elements ■ maximum norm ■ fully discrete • resolvent estimates • 
resolvent estimates • optimal error estimates • parabolic smoothing 


1 Introduction 

Let 17 be a Lipschitz domain in d = 2, 3 and I = (0, T). We consider the heat equation as a model of a 
parabolic second order partial differential equation, 

dtu{t,x) — Au{t,x) = f(t,x), {t,x) G I X f], 

u{t,x)=0, {t,x)£lxdn, (1) 

u{0,x) = uq{x), X G 

with a right-hand side / G L^{I\ LP{f2)) for some 1 < p, s < oo and uq G U'{Q), 1 < p < oo. 

The maximal parabolic regularity for uq = 0 says that there exists a constant C such that, 

\\dtu\\L‘(i-LP(0)) + \\Au\\L‘’(i-Lp{n)) <C\\f\\L-(i-Lp(0)), l<p,s<oo, for all / e L®(/; LP(17)), 

(see, e.g., llMT9ll20l ). The maximal parabolic regularity is an important analytical tool and has a number of ap¬ 
plications, especially to nonlinear problems and/or optimal control problems when sharp regularity results are 
required (cf. 021ll22ll23ll25]l l. Our aim in this paper is to establish similar maximal parabolic regularity results for 
time discrete discontinuous Galerkin solutions as well as for the fully discrete Galerkin approximations. Such 
results are very useful, for example, in fully discrete a priori error estimates and are essential in order to keep the 
spatial mesh size h and the time steps k independent of each other (cf. j^ ). In lIZTl we apply the results of this 
paper to establish pointwise best approximation estimates for fully discrete Galerkin solutions. 

Maximal parabolic regularity with applications to semidiscrete finite element Galerkin solutions in space were 
analyzed for smooth domains in II14II15II and for convex polyhedra in Il29l . Time discrete results are much less 
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known in the finite element community. Explicit methods are treated in Implicit Euler methods with 

pointwise norms in time are considered in 016II17I . A more systematic investigation of discrete maximal parabolic 
regularity for various time schemes was carried out by Sobolevskii and Ashyralyev and summarized in the book 

m. 

In this paper, we investigate maximal parabolic regularity for a family of time discontinuous Galerkin (dG) 
methods, which were first deeply analyzed for linear second order parabolic problems in ifTSll . There is a number 
of important properties that make the dG schemes attractive for temporal discretization of parabolic problems. Eor 
example, such schemes allow for a priori error estimates of optimal order with respect to discretization parame¬ 
ters, such as the size of time steps and the mesh width, as well as with respect to the regularity requirements for 
the solution (see, e.g., fllOlll 111 1. Different systematic approaches for a posteriori error estimation and adaptivity 
developed for finite element discretizations can be adapted for dG temporal discretization of parabolic equations, 
(see, e.g., Since the trial space allows for discontinuities at the time nodes, the use of different spatial 

discretizations for each time step can be directly incorporated into the discrete formulation, (see, e.g., Il3^ 1. Com¬ 
pared to the continuous Galerkin methods, dG schemes are not only A-stable but also strongly A-stable, (see, 
e.g., mi An efficient and easy to implement approach that avoids complex coefficients, which arise in the equa¬ 
tions obtained by a direct decoupling for high order dG schemes, was developed in llJTl . Eor the treatment of 
optimal control problems, Galerkin methods are particularly suitable since they expose an important property that 
the two approaches optimize-then-discretize, i.e., the discretization of the optimality system built up on the con¬ 
tinuous level, and discretize-then-optimize, i.e., discretization of the state equation and subsequent construction 
of the optimality system on the discrete level, lead to the same discretization scheme, (see, e.g., ||4]). Compared to 
continuous Petrov-Galerkin time-stepping schemes (see lIlSl for details), dG schemes also have the advantage that 
the adjoint state can use the same discretization as the state variable. This allows for unified numerical treatment 
and simplifies a priori and a posteriori error analysis, (see, e.g., imi32ll33]|34l ). 

The main results of this paper for the time semidiscrete discontinuous Galerkin Uk solution consist roughly of 
two parts. Eirst, for the homogeneous problem (i.e. / = 0) with uq G 1 < p < oo we show 




[Ttfcjm—1 


LP{0) 


c 

Cm, 


( 2 ) 


for m = 1, 2,..., M. Then, using this smoothing result, we also establish discrete maximal parabolic regularity 
for the inhomogeneous problem when uq — 0. We show, 


M 


M 


■,LP{n)) + \\^Uk\\L‘{I-LP{n)) + E 


\m—l 


\m—l 


['Cifcjm —1 


LP(f2)^ 


T 


(3) 

for 1 < s < oo and 1 < p < cxd, with obvious notation changes in the case of s = oo. In the case of the lowest 
order piecewise constant method, i.e., g = 0, the first terms on the left-hand side of the above estimates vanish. 
In contrast to the continuous case, the limiting cases s,p G {1, cxd} are allowed, which explains the logarithmic 
factor in Q. We also provide the fully discrete analog of (|2]i and Q. 

The rest of the paper is organized as follows. In the next section we introduce the discretization method and 
the resolvent estimates, which build the main analytical tool of the paper. For better communication of the ideas 
we first analyze the dG(0) method, which is technically much simpler, and in the following section we analyze the 
general dG(g) case. That is done in Sections |3]and|4l respectively. At the end of Section|4]we provide an example 
of how such maximal parabolic regularity results can rather easily lead to optimal order error estimates. Finally, 
Section |5] is devoted to fully discrete Galerkin solutions. In Section |6] we provide an extension of our results to 
the case of a general norm fulfilling a resolvent estimate. This generalization, being of an independent interest, 
is used, for example, in for derivation of pointwise interior (local) error estimates of fully discrete Galerkin 
solutions. 


2 Preliminaries 

To introduce the time discontinuous Galerkin discretization for the problem, we partition I = (0, T) into subin¬ 
tervals Im = of length km = tm - tm- 1 , where 0 = fo < G < • • • < tM-i < tM = T. The 
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maximal and minimal time steps are denoted by A: = max™ km and kmin = min^ km, respectively. We impose 
the following conditions on the time mesh (as in ED): 

(i) There are constants c,/3 > 0 independent on k such that 

^min ^ . 

(ii) There is a constant k > 0 independent on k such that for all m = 1, 2,..., M — 1 

1 km, 

K < < K. 


(iii) It holds k < jT. 

The semidiscrete space of piecewise polynomial functions in time is defined by 

Xl = {uk€L^iI;H^{n)): Uk\i^ €Vq{Hl{n)), m = 1, 2,..., M}, 

where Vq{V) is the space of polynomial functions of degree q in time with values in a Banach space V. We will 
employ the following notation for functions in 

M+ = lim u{tm + e), = liui u(tm-e), [u]^ = (4) 

£->■ 0 + £->- 0 + 

Next we define the following bilinear form 

M M 

B{u,(p) = ^{dtu,(p)i^xn + (Vu, V(p)/xr2 + '^{[u]m-i,‘Pm-i)n + (uo ,¥^ 0 )^, (5) 

m—1 m—2 

where and xt? are the usual space and space-time inner-products, is the duality product 

between L'^{Im', H~^{^)) and L^{Im', note, that the first sum vanishes for u G The dG(g) 

semidiscrete (in time) approximation Uk G of dB is defined as 


B{uk,(pk) = {f,y:’k)ixn + {uo,(p'l^_Q)n for all tpfc G X®. (6) 

Rearranging the terms in Q, we obtain an equivalent (dual) expression of B: 

M M-1 

B{u,(p) = - ^ {u,dtip)i^xn + (Vu, Vvj)/xr 2 - ^ (u“, [‘p]m)n + {uM,‘PM)n- (7) 

m—1 m—1 


The analysis of such schemes in non-Hilbertian setting is usually done by using a semigroup approach that 
represents time stepping formulas as a Dunford-Taylor integral in the complex plane Bdl Ch. 9]. This approach 
requires certain resolvent estimates. For Lipschitz domains and a given 7 G (0, 7 r/ 2 ), the resolvent estimate (see 
HOl l guarantees the existence of a constant C such that for all u G L'p{Q), 1 < p < 00 , and any z G C \ Ti’.y the 
following estimate holds: 

||(z-l-Z\) ^M||LP(r 2 ) < Y^p-|^||u||LP(r 2 ), ( 8 ) 

where the Laplace operator — Zi is supplemented with homogeneous Dirichlet boundary conditions, and 

= {z G C : |arg(z)| < 7 }. (9) 

Using the identity A{z A)~^ = Id —z(z -I- A)~^, one immediately obtains, 

\\A{z +A)~^u\\LP{a) <C\\u\\lpi^q), zgC\X^, 1 < p < 00 , itGL^(l7). (10) 

We note, that all our results for semidiscrete solutions hold if we replace the Laplace operator — A with a more 
general self-adjoint second order elliptic operator A provided it satisfies (O. 
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3 Estimates for dG(0) 

For the ease of the presentation, we first establish the results for the lowest order piecewise constant discretization 
dG(0). In this case, we use the following notation, 

Uk,m=Uk\l^, ul.^=Uk,m+l, 'U‘k,m=Uk,m, 771 = 1, 2, . . . , M - 1. (11) 

First, we establish results for the homogeneous problem. In this case the dG(0) method is equivalent to the Back¬ 
ward Euler method. 


3.1 Results for the homogeneous problem 


Let f = 0, uq € LP{f2) and let Uk S be the semidiscrete approximation of ([T]) defined by 

Biuk,Xk) = {uo,Xk.,i), Vxfe G (12) 

i.e., the dG(0) solution Uk satisfies 

Uk,i — kiAuk,i = uq, 

'^k,m kraAuk^m — '^k,m—li 777 = 2, 3, ... , iff. 

The first result shows that the solution can not grow from one time step to the next one. 

Lemma 1 Let Uk be the solution of (I12i . Then, for uq G 1 < p < oo there holds 

\\uk,m\\LP{0) < ||?^o||LP(t2) V?77 = 1, 2, . . . , M. 

Proof First, we assume uq G and establish 

l|Mfc.m||L~(r7) < ||uo||l~(j7) 777 = 1, 2, . . . , M. (14) 

It is sufficient to consider only a single time step. 


Uk,i-kiAuk,i=uo. (15) 

We want to show that < 11 770 Assume it is false. Let xq G 17 be a point where attains a 

maximum. By ifTSl Theorem 3.3], we know that Uk,i G (70(17), hence, there exists an open ball Bs{xq) of radius 
(5 > 0 centered at xq with Bs {xq ) C 17 such that 


Uk,i(x) > ||wo||L~(t 2 ) for all x G Bs(xo). 


Hence, 

By the maximum principle, from 


Uk,i(x) - uo(x) > 0 onBs(xo)- 


-Auk,i = — (uo - Uk,i) < 0 onBs(xo), 
ki 

we obtain a contradiction to the assumption that Uk,i has a maximum at the interior point Xq- This contradiction 
establishes (Eli. Next, using a duality argument, we will show 


||ufc,l||Li(t2) < ||Mo||Li(t3)- (16) 

Consider the problem, to find Zk.i G ffp (17) that satisfies, 

Zk,i - kiAzk,i = zo, with zo = sgnuk,!- 


The solution Zk.i can be thought of as a single step of the dG(0) method to a parabolic problem with initial 
condition sgnufc,i. Thus, 

l|wfe,l||Li(t3) = (uk,l,Zo) = (zk,l,Uk,l) + ki(Vzk,l,VUk,l) = (uo,Zk,l) < ||Mo||Li(l7)ll^o||L“(r2) < II^o||l1(17), 

where we have used (fT4l i for Zk and the fact that || 2 o||l°°(i 7 ) = llsguTtfe^i= 1. This establishes (fTfil) . 
Interpolating, we obtain the lemma for 1 < p < c». 
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Next we will establish a smoothing result. 

Theorem 1 (Homogeneous smoothing estimate) Let Uk G be the solution of (fTST i with uq g LP{f2), 1 < 
p < oo. Then there exists a constant C independent ofk such that 

C 

\\^Uk,m\\LP(Q) < — ||?^o||LP(t2)> m = 1,2, . . . , M. 

Proof The proof is given on page 1321 in Ha for the norm, but the proof is valid for the L'p{Q) norm as 

well by using the resolvent estimate ® with respect to the L'P{f2) norm. 

Remark 1 Let Uk G be the solution of (fT^ with uq G LP{S1), 1 < p < oo. Then there exists a constant C 
independent of k such that 

||'ttA;,m||LP{t7) “b {im ^/) || || ||LP{t7) ; V?7t > / ^ 1. 

From (foi l, we immediately obtain the following result. 

Corollary 1 Let uu G be the solution of (112b with uq G L^(f2), 1 < p < oo. Then there exists a constant C 
independent ofk such that 


[ttfcjm — 1 


LP{n) 


C 

< ^l|MolUp(t2), m 

Lrn 


1,2,...,M. 


3.2 Results for the inhomogeneous problem 

Now we consider Uk G X^ to be the dG(0) solution to the parabolic equation with uq = 0, i.e., Uk satishes, 

B{uk,Pk) = {f,Pk)ixn, 'ipkGXl. (17) 

Thus, the dG(0) solution satishes 

Uk,i - kiAukp = kifi, 

(lo) 

Uk,m = 2, 3, . . . , TV/, 

where 

/m(-) = / f{tr)dt. 

JI j.ri 

Since fm is the Lf projection of / onto the piecewise constant functions on each subinterval Im, we have 

max \\fm\\LP{n) < C\\f\\L^(i.^LP{n)), 1 < P < oo, 

M 

XI l<P<oo, l<r<00 

m—1 

We now state our main result for the dG(0) approximations. 

Theorem 2 (Maximal parabolic regularity) Let 1 < s,p < oo and uq = 0. Then, there exists a constant C 
independent ofk such that for every f G L^{I; and Uk satisfying (117b . the following estimate holds: 

T 

||^Mfe||L-(/;LP(n)) < C'ln— 1 <S<CX), l<p<00. 


(19a) 

(19b) 
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Proof Using Cl, we can write the dG(0) solution as 

m / m — /+1 \ 

'^k,m — ^ ^ I j + 1 ^) I fh ^ — I? ^5 • ■ ■ 1 ^ 

/=! \ 1 = 1 / 


where r(z) = (1 + z) Then, 

m / m—Z+1 \ 

= Z\ r(-A;^_j+iZ\) /;, m = 1,2,..., M. 

7=1 \ 1=1 / 


Hence 


m 

||Z\Mfe ||lp(J7) ^ ^7 


7=1 


m —7+1 

^ n ?’(-+-j++) j fi 


7 = 1 


m 




\ ' / LP{0) 

From Remark[T] since each term in the sum on the right-hand side can be thought of as a homogeneous solution 
with initial condition fi at t = f;_i, we have 


m—7+1 

^ n r{-km-j+i^) j fi 


< 


C 


Lp(n) 


^1 — 1 


\\fi\\LP{n)- 


Thus, we obtain 


rn J 

|+Mfc,m||LP(t2) < C* ^ - -—^7- \\fl\\LP(0), 171 = 1,2, ... ,M. 

^7 — 1 


where in the last step we used that 


tm — tl-l 


ptm-l 


< 1 + 


dt 


= l + ln^<Cln^, 


'0 i 


by using the assumption +iin > Ck^ and k < ^. 
For s = 1, we have 


M 


M m 

\\^Uk\\mi-,LP(n)) = E km E km Er^ 77 — \\fi\\LP(o)- 

m—1 m—1 1—1 ^ ^ ^ 

Changing the order of summation and using ( fT9] t, we obtain 


M 


M ivi J 

\\^uk\\Li(i-,LP{o)) < c'y]fciii/iiiip(j7) ^ 

, - , 


7=1 


n—l 


U-1 


M 


rp rp 

- E ki\\fi\\LPiO) E C'ln — 

^ 7=1 


where we used again that 


M 


n—l 


T 

■ < Cln-. 
^7 — 1 k 


i:+ 


Interpolating between s = 1 and s = oo, we obtain the result for any 1 < s < oo. 


( 20 ) 


1 — 1 

For s = oo, we obtain from the above estimate and using d, 

m 

|+Mfe||Loo(j.ip(^2)) = max \\Auk,m\\LPin) max - - — WfiUpiO) 

l<m<M l<m<M “—^ tm — ^7 — 1 

^ k T 

< C^m|x7l/,||„,„, +^77 S C'>"t ll/ll win)). 


( 21 ) 
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Remark 2 The appearance of the logarithmic term is natural for the critical values s = 1, p = 1, s = cxd, or 
p = oo, since the corresponding maximal parabolic regularity results for the continuous problem hold only for 
1 < s,p < oo. For s = 2 or p = 2, the power of the logarithm can be lowered. Thus, for p = 2, from ll^ we 
know, 

and from (l20l i. we have 

T 

II^Ufc||L“(/;L 2 (r 2 )) < Cln —1 < S < CX). 

Interpolating between s = 2 and s = oo and between s = 2 and s = 1, we obtain 

l«- 2 | 

3 

l|/l|L»(/;L2(r2)), for any 1 < s < OO. 


\\^Uk\\L‘(i-L'^(n)) < C" ( In-^ 


T 


Similarly, we can obtain. 


< C I In 


IP-2| 


\\f\\L^i-LP(n)), for any 1 < p < oo. 


Corollary! (Maximal parabolic regularity for jumps) Let 1 < s,p < oo and uq = 0. Then, there exists a 
constant C independent of k such that for every f S L‘^(I; LP[L2)) and Uk satisfying (I171 l. the following estimate 
holds, 


max 

l<m<M 




LP(i7) 


< C'ln^||/||icx,(j.ip(j 7 )), l<p<oo, 


M 


I] fcr, 


[Uk] 


m—1 






‘iI;LP(a))j I < S < OO, l<p<oo, 


where the jump term [rtfejo at t = 0 is defined as Ukp- 
Proof Since by (fTsT i on each time subinterval Im we have 

—1 — ^tlk,m “f fm, Tfl = 1,2,..., AT, 

by using Theorem|2l we have 


T. 


max ^ \\[uk]m-i\\LPm) < raax {\\Auk,m\\LPin) + ll/m||LP(r3)) < C'ln-p||/||L=o(7;LP(r2)). 

Cm <CM ^ > 1 <'m.< A// K* 


l<m<M ^ l<m<M 

Similarly, using Theorem]!] for 1 < s < oo we have 

M r 1 S M 

[Uk\m-1 


m—1 


LP{f2) m—1 


Auk.mWlpf^Q) + ll/mll^p^o^ ) < Cg ( In 


\LPiO) ^ 


IL'>(7;Lp(J7))i 


where the constant Cs depends on s. By taking the s-root we obtain the corollary. 


4 Estimates for dG(q) 


In this section we will establish the dG(q) version of the results from the previous section. It is convenient to 
introduce some additional notation. Let q > 1 and ipiit) G F’g([0,1]), I = 0,1,..., g be the standard Lagrange 

basis functions on the interval [0,1], i.e., where Sij is the Kronecker symbol. Then for any Uk G 

on the time interval Im = (fm-i, fm] we have 


i =0 


(22) 


with (7™ G Uq (17) independent of t. In this notation, we have 


= and ul^ = U^. 



















8 


Dmitriy Leykekhman, Boris Vexler 


4.1 Results for the homogeneous problem 


Let Uk G X^. be the semidiscrete in time solution to the parabolic equation with / = 0, i.e., 

B{uk,ipk) = V(/7fe G X’. (23) 

Alternatively, on a single interval Im, we have 

Ul = ri^o{-kiA)uo, l = 

Ur = ri,oi-kmA)U^-\ Z = 0,l,...,g, m = 2,3,...,M, 
where the rational functions r; o are of the form, 

n,o(A) = ^^, l = 0,l,...,q, (25) 

P(A) 

with p being a polynomial of degree q + 1 with no roots on the right-half complex plane and pi^, I = 0,1 ,..., g 
being polynomials of degree q (cf. IfT^ . page 1322). Since o(A) is a subdiagonal Fade approximation of e“^, 
we also have (cf. ii) 

?'9,o(0) = Pq, oiO) = p(0) = 1 and |rg,o(A) - e“^| = 0(|Ap«+^), (26) 


as A —>• 0. The rational functions r; g satisfy the following properties, which we will often use 

n.o(O) = 1, and n,o(A) - 1 = / = 0, l,...,g, (27) 

P(A) 

where pi{X) are some polynomials of degree q. The first property follows, for example, by considering the homo¬ 
geneous Neumann problem with initial condition ug = 1. Then the exact solution u and the dG{q) solution Uk 
are the same and equal to 1, i.e., u = Uk — I- Hence, all nodal values [//" = 1 for all m = 1,2,..., M and 
Z = 0,1,..., q. For example for m = 1, we have 

l = Ul = rip{-kiA)uo = rip{-kiA)l = ri,o(0), 

and as a result r;_o(0) = 1. The second property in (l27l) is just a consequence of the first one. 


Remark 3 The dG(l) solution Uk on each subinterval Im is of the form 



and the rational functions are p(A) = 1 -f |A -f rg o(A) = 1-1- |A, and ri g(A) = 1 — -1. 


For later proof we require two supplementary results. 

Lemma 2 Let the rational function r{z) be of the form r(z) = where p{z) is a polynomial of degree q -\- 1 
with no roots on the right half complex plane and p{z) is a polynomial of degree q, for some q > 0. Then, there 
exists a constant C independent ofk>0, such that for any g G 


r{-kA)g\\Lvi^n) < C\\g\\Lv(n). 


(28) 


Proof For simplicity we assume that the roots zi, Z2-, ■ ■ ■, Zq of p are pairwise distinct. If it is not the case, the 
argument can be slightly modihed. For g = 0 we have r{z) = 73 ^ and the desired estimate follows directly by 
the resolvent estimate ®, since 

'■(-“)’ =-|(t + ^) 

and therefore by ® 
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For 5 > 0 we use the partial fraction decomposition 


r{z) = 


Ci 


q 

V — 

z — z, 

i=0 


with some Ci e C. Applying the estimate for to each summand we obtain 

Wghno), 


which completes the proof. 

Lemma 3 Let the rational function r(z) be of the form r{z) = where p{z) is a polynomial of degree q+1 

with no roots on the right-half complex plane and p{z) is a polynomial of degree q, for some g > 0. Then for any 
g G L'P{f2) with Ag G L'p{L 2), \ < p < oo, there exists a constant C independent of k such that 


\\r[-kA)g\\LP^Q) < Ck\\Ag\\LP^^Q). 

Proof This lemma is just a consequence of the previous one. We set f(z) = and obtain: 

r{—kA)g = —kAf{—kA)g = —kf{—kA)Ag. 

The the result follows by Lemma|2] 

Lemma 4 Let the rational function r{z) be of the form r{z) = where p{z) is a polynomial of degree q + 1 

with no roots on the right half complex plane and p{z) is a polynomial of degree q, for some q > Then, there 
exists a constant C independent of k, such that for any g G LP{f2) 


\\r{-kA)g\\LP(o) <C\\g\\Lp{n)- (29) 

Proof We set f(z) = and obtain: 

||r(-fcZ\)g||ip(f2) < k\\Af{-kA)g\\Lp(^Qy 

The estimate 

\\Af{-kA)g\\LP(^a) < j\\g\\LP(0) 

is provided on the top of page 1322 in ifT^ using a decomposition r(z) = ri(z) + r 2 (z), where ri(z) = 
with zq being a root of p(z) and c such that the degree of the polynomial in the numerator of r 2 (z) is less or 
equal q — 1. Then the estimate for Ari{—kA)g follows directly by applying a dG(0) type argument and the term 
Ar 2 {—kA)g is estimated using the Dunford-Taylor formula. 


Next we provide some properties of the dG(q) solutions of the homogeneous problem. 

Lemma 5 Let Uk be the solution of (1231) with uq G L'P[f2), \ < p < oo. Then, 

ll'tifc||L°°(/„;LP(l7)) < C'||uo||lp(i7), Vto = 1, 2, . . . , M. 

Proof The proof is given in ifT^ Thm. 5.1] for the Lf{Q) norm, but the proof is valid for the LP(17) norm as well 
by using the resolvent estimate ([8]) with respect to the L^’(l7) norm. 

Theorem 3 (Homogeneous smoothing estimate) Let Uk be the solution of (|2^ with uq G Lp{L2), 1 < p < oo. 
Then there exists a constant C independent of k such that 

C 

< — ||uo||LP(r2), m = 1,, 2 ..., M. 

Proof Again the proof is given in ifT^ Thm. 5.1] for the Lf{Q) norm, but the proof is valid for the LP{L2) norm 
as well by using the resolvent estimate (|^ with respect to the LP{f2) norm. 
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Remark 4 Notice that the statement of Theorem[3]is equivalent to 

\\AUr\\Lv(a) < ^\\uo\\L^(n), m = Z = 0,l,...,g, (30) 

which we will use in the following proofs. 

Remark 5 Let Uk be the solution of (l23t . Then there exists a constant C independent of k such that 

\\ulJ\\LP{,n) + {tm - tn)\\Auk,m\\L--(Im-,LP(n)) < II Lp (t?), m > Tl, n=l,2,...,M, 

or in terms of nodal values 

\\U^\\LPin) + {tm-tn)\\AUr\\LPin)<C\\U^hpin), m > n, n = l,2,...,M, I = 0,1,... ,q. (31) 

Theorem 4 (Homogeneous smoothing estimate for jumps) Let Uk be the solution of (|2^ with uq € LP(f7), 
1 < p < oo. Then there exists a constant C independent ofk such that 

C 

< —l|Mo||LP(t 2 ), m = 

LP{n) 

where [Mfe]o = Uq - Uq. 


[ttfcjm —1 


Proof For m> 1, using dH, we have 


Kl^-i = = rQ,Q{-kraA)U'U-^ - = M-k^A) - ld)U^-\ 


Using dZTl) and Lemma[3l we obtain 


[tifcjm —1 


LP(a) 




Now by Remark|4]and the assumption on the time mesh {ii), we obtain 

- 7-lko||LP(t2) < 7—ll^‘o||LP(t2)- 

' ^ ^m—1 


That finishes the proof for this case. 
For m = 1, by Lemma|5]we have, 


[Mfc]o 

ki 


Lp{n) 


— ||(7o - tio||LP(r2) < ■^||wo||Lp(r2) = — ||mo||lp(j7)- 


Similarly, we can obtain the corresponding result for the time derivative. 

Theorem 5 (Homogeneous smoothing estimate for time derivatives) Let Uk be the solution of (|2^ with uq G 
L'P[f2), 1 < p < 00 . Then there exists a constant C independent of k such that 


C 

\\dtUk\\L°-{I^-LP{0)) < —|| mo || lp ( 17 )- 

Lm. 


Proof For m > 1, using ^ and da, we have 


dtUk\i^ = k^ ^ Ur{x)'tlj'i 


1=0 


t tm — 1 


= k^^nfi{-kmA)%l;'i 


1=0 


t im—l 


U^-^x). 


By the fact that I'" ^ ) = 1 we have ''Pi = 0. Using dZTl l. i.e., r; o(0) = 1 we obtain 



zpt{z) 

p{z) 
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where p{z) is the same polynomial as in (l25l l and pt (z) is some polynomial of degree q — 1 whose coefficients are 
time dependent, but uniformly bounded on 1^- Thus again by Lemma [3 we obtain 

\\dtUk\\L'^{i^-LP{n)) < C\\AU^~^\\LP{n)- 

Remark]?] and the assumption on the time mesh (ii), finishes the proof for to > 1. 

For TO = 1, by Lemma]5]we have, 

^ ^ (j 

\\dtUk\\L°°(ir,Lp{n)) < Ck^^'^\\Ul\\LP(n)Wi\\L°-{h) < — || mo || l !>( j 7 )- 

1=0 ^ 


4.2 Results for the inhomogeneous problem 


In this section we establish properties of the dG(g) solution Uk G to the inhomogeneous parabolic equation 
with uq = 0, that satisfies, 

B{uk,(pk) = {f,^k), 'ipk&Xl. (32) 

Alternatively, on a single time interval Im, we have 

(33) 

Ur = np{-k^A)U^-^ + km^ri^,i-kmA)f^, Z = 0,1,... ,q, TO = 2,3,... ,M, 

3=0 


where 


and the rational functions 


/r(-) 



fit, ■)'<P3 


t fm —] 


dt, 


n,3 = 


pijW 

pW 


l,j = 0, l,...,g. 


(34) 


are as in the homogenous case with p being a polynomial of degree g + 1 with no roots on the right half complex 
plane and pij, l,j — 0,1, ■ ■ ■, q being polynomials of degree q (cf. Ill2l,page 1322). 

Notice that for to = 1,2,..., M, 


ll/]"llLP(t2) < C'll/llL“(/m;LP(J7)) ^nd \\f^\\LP(n)<Ckm^\\f\\L^I„,-,LP(a))- (35) 

Theorem 6 (Maximal parabolic regularity) Let Uk satisfy (l32l i with f € L®(/; LP{Q))for 1 < s,p < oo. There 
exists a constant C independent ofk and / such that 

T 

\\Auk\\Lp{i;LP{n)) < C!ln — \\f\\k‘,{i-LP(n))- 

Proof Using (l33l l. we have the following representation 

m—1 j m—n—1 \ 

ur = kmGT + ryo{-kmA) ^ fcJ n U^oi-km-j-iA) (36) 

n=l \ j=l I 


where 




J2n,r-kmA)fr, 


3=0 


TO = 1,2,..., M. 
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with the usual convention that 0^=1 ™ empty product. The proof now follows along the lines of Theorem|2l 

Taking the Laplacian of both sides we obtain 


m —1 / m—n — 1 


Aur = kmAGT + An^^{-kmA) ^ H r,^o{-km-j-iA) | G” 

n=l \ j=l 


and as a result 

WAUrU^^n) < Wk^AGTWL^^n) 

By Lemma|4] we have 


m —1 / m—n — 1 

Ari^o{-kjnA) ^ j rq^o{-kni-j-iA) ] G^ 

n=l \ j=l 


Lp{n) 


\\kmAGT\\LP(n) < C max Wf^hpin), / = 0,l...,q, 


and by Lemma|2]we also have 


l|GrilLP(r?) <C maxJ|/™|U.(^), l = 0,l...,q. 


(37a) 


(37b) 


On the other hand by Remark |5] for any I = 0,1,..., g, since each term in the sum on the right-hand side can be 
thought of as a homogeneous solution with initial condition G^ at f = <„_i, we have 


m—1 / m—n—1 

Ari^o{-kmA) ^ I rq^o{-km-j-iA) | G^ 


n—1 


To establish the result for s = oo, we observe 


m—1 


k 


<G^ ||G"|Up(r,). (38) 

^m ^n — 1 


LP{0) 


n—1 


\\^Uk\\LP-iI-LP{n)) = 


<C max max Wf^WiPio^-^C max > 

l<m<M0<j<g ^ ^ 1 <rrr><'/\/f 


l<m<Af tjji tji—i 

n—1 


\G 


q \\LP(n) 


<C max max ||/’"||ip(f 2 ) 1-f max 

l<m<M 0<j<q ^ ^ M T<^<M 


l<m<M ^ t^n tn—1 
n—1 y 


T 


<Gln-- max max \\frW lp( o), 

where in the last step we used (EB. Using drsl l we can conclude that for s = oo 

T T 

\\Auk\\Lp^(i-LP{n)) < < C\n—\\f\\L,^(i.LP(a)). 


Similarly, for s = 1, we have 


M 


\\-^Uk\\L^i-Lp{n)) < ^ maxJ|Z\C/n|iP(r 3 ) 


m—1 

M 


M m — 1 


<GY,k^ max Wfphp^a) + E E 


m—1 

M 


m—1 n—1 

M m — 1 


< G E max \\fp\\LP{n) + C* E E 


m—1 

M 


m—1 n—1 


^n — 1 

kn 

^n — 1 ^ 


max ||/"||LP(f 2 ) 


m=l n=l 
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Changing the order of summation and using (1211 1 we obtain, 

Mm-, M 


M 




0<j<9 II II • 


M 


Thus, by using (1351 >. we have 


M 


T . . T 

\\^Uk\\L^{i-LP(n)) <C\n—'^km \\fJ^\\LP{n) < C'ln —||/||ii( 7 .ip(j 7 )). 


Interpolating between s = 1 and s = oo we obtain the result for any 1 < s < oo. 

Remark 6 As in the case of dG(0) the appearance of a logarithmic term is natural, since in contrast to the contin¬ 
uous case the choices s, p G {1, cjo} are allowed. The power of the logarithm can be improved for p = 2 or s = 2. 
In fact, we can obtain the following estimates. 


||^W/c||l'*( 7;L2(J7)) < C* I In 


°- 2 | 


and 


|L=(/;L 2 (n)), 


II^Wfc||L2(/;iP(72)) < C I In — I ||/||l2(/;Lp(17))- 


, lp-^l 

T\~l^ 


Theorem? (Maximal parabolic regularity for jumps) Let Uk satisfy (l32l) with / G L’^{I\Lp{Q)) for 1 < 
s,p < oo. Then there exists a constant C independent of k and f such that 


max 

l<m<M 




M 




h 

[tifejm —1 


T 


Lp{n) 


< C'ln--||/|| 7 ,oo( 7 . 7 ,p( 72 )), fors = oo, 


LP{Q). 


T 

< Cln—|l/||i»(/.LP( 73 )), forl<s<oo. 


Proof Using (l33l l and (l36l) . we have the following representation for the jump terms 

[ukU-i _ 

k h 

= G™ + {ro,o{-kmA)U^-^ - =G^ + k-^ (ro,o(-fc™Zl) - Id) U^-\ 

Using that rop — 1 satisfies (l27l i and using Lemma[3l Lemma|2] and proceeding similarly to the proof of Theorem 
1^ we have 

kf^\\[uk\m-l\\LP{n) < G (||G™||iP(72) + ||Z\C/™“^||iP(j7)) 


m—1 




^n —1 


max Wf^hpio) 


( 39 ) 


- ‘^'12 7 - 7 - \\f^\\Lp{n)- 

, tm — tn-1 0<3<q 
n—1 

Now, the proof of the cases 5 = 1 and s = oo is identical to the one of the previous Theorem[6]and we have 


max 

l<m<M 


M 

m—1 




h 

^m 


LP{0) 


LP{f2) 


< C'ln-||/||i=o( 7 .pp( 77 )), l<p<cx), 

T 

< Gln--||/||pi( 7 .pp( 72 )), l<p<oo. 
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For 1 < s < oo using the Holder inequality with - + ^ = 1, we obtain, 


—1 ^ II mil 

< C' :- max \\fn\LP{0) 

LP{n) ^n-l 0<j<q 


1/s 




l/s 


^n—l 


im tn—l 


(40) 


< C ( In ? 


I/s' 


E 


l/s 


^n — l 


max Wf^Wh/a) 


Hence 


M 


E 

m=l 




< C In 


LP{n) 


y \ s/s' M 




m=l n=l 


Changing the order of summation, we obtain 


M 


E 




LP(r 2 ) 


rji\ s/s' M 
n—1 

1+s/s' M 


< C In 


0<]<q ^ tm- tn-1 


M 

E 


E^"n’^f<^ 11^: 


n=l 


0</<9 


= C7 In 


T 


\L‘{i-,Lp(n))- 


Taking the s-root we finish the proof. 

Theorem 8 Let Uk satisfy (132b . Then there exists a constant C independent of k and f such that 


M 


I^E ^ C'ln-||/||i.(/.iP(f 2 )), l<s<oo, l<p<oo. 

Proof Similarly to the proof of Theorem|4] using (l22li and (l33b . we have 

dtnk\i^=kf,^j^Ur{x)iP[ +j^GT{xW, 

7 _n \ f^m / \ f^TTL / 


Z=0 

q 


= U^-\x) + j2^T{xWi 

1=0 \ m / 


/ / f ^m — 1 


Using dZTl) and X]?=o V'; ~ conclude that 

, fzpt{z) 


^Xko{z)'f'i 


1=0 


p{z) 


where p{z) is the same polynomial as in (l25l l and pt{z) is some polynomial of degree q whose coefficients are 
time dependent, but uniformly bounded on Im- Thus again by Lemma|3]and Lemma|4] we obtain 

\\dtUk\\L'^{i^-,LP{a)) < C\\AU^~^\\LP{n) + C^max^ ||/J”||LP(n)- 


The rest of the proof is identical to the proof of the previous theorem. 
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4.3 Application to optimal order error estimates. 


As an application of the maximal parabolic regularity, we show optimal convergence rates for the dG(( 7 ) solution. 
First, we establish that the error is bounded by a certain projection error. A similar result was obtained for the 
L^(/; L^(J7)) norm in ll^ . First, we define a projection tt^ for u G C'(/, L^(J7)) with Trku\i^ G Pq{L‘^{n)) for 
771= l,2 ,...,Mon each subinterval Im by 

(TTfeU - 77, (/))/,„ X 72 = 0, Mcj) € Pq-l{Im,L^{^)), 9 > 0, (41a) 

(41b) 

In the case 9 = 0, Tr^Tt is defined solely by the second condition. 

Theorem 9 Let u be the solution to ([T]i with u G C{I', LP{f2)) and Uk be its dG(q) approximation for 9 > 0. 
Then there exists a constant C independent ofk such that 

T 

\\u - Uk\\L-iI-LP(n)) < C'ln — 1|77 - 7rfc77||L«(/.Lp(f2)), 1 < s,p < CX), 

where the projection tt^ is defined above in (141b . 


Proof Put e := u — Uk = {u — Tr^it) + (Tr^Tt — Uk) '■= rjk + Cfc- For 1 < s,p < c», we have 


\^\\l‘(I-,Lp(0)) 


sup 

(n)) 

ll^ll / / = 

' "ls (17)) 


(e,V') 


7x72) 




For each such tl>, we consider a dual problem for Zk G satisfying 

B{pk,Zk) = iPk,i’)ixO forall G 


1 

P 



Thus, we have 


(e,7A)7x72 = {'nk,f)ixn + {^k,'f)ixn ■= Ji + ^2- 

Using the Holder inequality, we find 

Jl < l|t?fe||L»(7;LP(72))ll'0llL'»'(7;7!>'(72)) ^ Ihfc ||l'= (7;Lp(72)) • 

On the other hand using that B(u — Uk,Xk) = 0 for any Xk G X^, we obtain 

M 

J 2 = B{^k,Zk) = -B{rik,Zk) = ^ {'nk,dtZk)i^xn - {'^r]k,'^Zk)i^xn + [2:fe]m)72 

m—1 

= —{Xrjk, Xzk)ixn, 

where we used that the first sum vanishes due to (14lab and the sum involving jumps due to (I41bb . Integrating by 
parts in space, using the Holder inequality and Theorem we obtain 

J 2 = -i'^rik,'^Zk)lxO = iVk, 2i\Zk)lxQ < ll9fc||L'’(7;L!>(72))l|4\^fc||is'(7.iP'(72)) 

T T 

~ ll^'=llr,“(7;LP(72))ll'0llL'>'(7;L!>'(72)) ^ C'In||??fe || L«(7 ;Lp(72)) • 


Combining the estimates for Ji and J 2 we obtain the result. 

If the exact solution is sufficiently smooth then the above result easily leads to an optimal convergence rate, 
modulo a logarithmic term. 

Corollary 3 Let u G (/; LP(f2)) be the solution to dTJ and Uk be its dG(q) approximation for 9 > 0. Then 

there exists a constant C independent ofk such that 


T 

u — Uk\\L‘'{I-,LP(a)) < Ck'^^^ In — ||77||vy,+l,«(7.ip(f2)), 1 < S,p < CXD. 
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Remark 7 The above result can be extended to the case of non-homogeneous Dirichlet boundary conditions. Let 
g € C{I; L‘^{Q)) n L^(/; i?^(l7)) be given and consider the equation 

dtu{t,x) — Au{t,x) = f(t,x), {t,x) € I X fl, 

u{t,x) = g{t,x), (t.x) G I X ^^2, 

u(0,a;) = uo(x), x G f2. 

It turns out, that it is convenient to use TT^g as boundary conditions for the semidiscrete solution, i.e. 

UkGnkg + Xl : S(ufc, </5fc) = (/, <Pfc)/xr2 + (mq, for all G 

Then following the lines of the proof of Theorem|9]and using that — Uk has homogeneous boundary 

conditions, i.e., ^k G X^, we obtain 


(Cfe,'</’)/xi7 =-(Vpfe, Vzfe) = (77fe,Z\zfc)/xr2 + / / {g - Trkg)dnZkdsdt. 

Jl JdO 

Under an additional assumption on 17 that for any v G iTg (17) with Av G V (17) the estimate 

ll^n^llLp'(9r2) — 

holds, we obtain 

T 

||u - Uk\\L‘(I-LP{n)) <Ch\ — (||m - ■Kku\\L‘{I-LP(Q)) + \\g - '^k9\\L‘‘{I;LP{dO))) , l<S,p<00. 

The above assumption is fulfilled, for example, if on 17 the P elliptic regularity holds. 


5 Fully discrete solutions 

In this section, we consider the fully discrete approximation of the equation ([T]). From now on we assume that 
the domain 17 is a polygonal/polyhedral convex domain. For h G (0, /ig]; > 0, let T denote a quasi-uniform 

triangulation of 17 with mesh size h, i.e., T = {t} is a partition of 17 into cells (triangles or tetrahedrons) r of 
diameter hj- such that for h = maxr hr, 

diam(T) << CItI^, Vr G T, d=2,3, 

hold. Let Vh be the set of all functions in that are polynomials of degree r on each r, i.e., 14 is the usual 

space of conforming finite elements. To obtain the fully discrete approximation we consider the space-time finite 
element space 

^1% = '■ Vkh\im ^’PqiVh), rn= 1,2,.. .,M, g > 0, r > 1}. (42) 

We define a fully discrete analog Ukh G X^’J^ of Uk introduced in (| 6 ]l by 

B{ukh,Pkh) = if,Pkh)ixn + {uo,(Pkh)^ for all (pkh G X^’l. (43) 

Moreover, we introduce the discrete Laplace operator Ah ■ Vh ^ Vh by 

{-AhVh, x)n = (Xvh,Xx)a, Vx G 14. 

The semidiscrete results from the first part of the paper translate almost immediately to the fully discrete setting 
provided we have the corresponding resolvent estimate, 

\\{z + Ah)~^x\\Lp( 0 ) < ^ ||xIIlp(j?), VzGC\r.^, MxGVh, l<p<oo, (44) 

with some constant C independent of h. Such a result was established in 13 for smooth domains. Later it was 
extended to convex polyhedral domains in ll^ (for some 7 > 0 ) via stability and smoothing properties of the 
semigroup Eh{t) = and directly for an arbitrary 7 > 0 but with logarithmic dependence of the constant C 

on h in ll26ll . 




Discrete maximal parabolic regularity for Galerkin finite element methods 


17 


5.1 Result for the homogeneous problem 

Let Ukh G be the fully discrete dG((?)cG(r) solution to the parabolic equation with / = 0, i.e. 

B{ukh,‘Pkh) = {uo,(pth,o)’ (45) 

Theorem 10 (Fully discrete homogeneous smoothing estimate) Let Ukh be a solution of (l45l l with ug € Lp({ 2), 
1 < p < oo. Then there exists a constant C independent of k and h such that 

C 

\\dtUkh\\L°^(I^-LP{Q)) + \\^hUkh\\L’=°(I^-LP(Q)) + \\[ukh]m-l\\LP{Q) < — 11^0 ||L!>(fi) > 

for TO = 1, 2,, M. 

5.2 Results for the inhomogeneous problem 

Let Ukh S ^k'^i dG(( 7 )cG(r) solution to the inhomogeneous parabolic equation with uq = 0, i.e. 

Bi^UkhiT^kh) — tf^Pkh^j '^Pkh ^ ^k,h' (^^) 


Theorem 11 (Fully discrete maximal parabolic regularity) Let Ukh satisfy (l46l i with f G L^{r,LP{fI)), 1 < 
s,p < oo. Then there exists a constant C independent of k and h such that 

1 1 


M 


M 


^ ^ \\^tUkh\\h‘{I,y^-,LP(Q)) 1 ~^\\^hUkh\\L‘{I-,LP(Q))~^\ ^ ^ 


[ti/c/i] 


m— 1 


T 

< Cln-I 


LP{n)^ 


lL = (/;LP(t3))) 


with obvious notation changes in the case of s = oo. 


5.3 Application to optimal order error estimates. 

Similarly to the semidiscrete case, as an application of the maximal parabolic regularity, we show optimal conver¬ 
gence rates for the dG(g)cG(r) solution. 

Theorem 12 Let u be the solution to ([T]) with u G C{I; L^(J7)) and Ukh be the dG(q)cG(r) solution for q > 0 
and r > 1. Then there exists a constant C independent of k and h such that for 1 < s,p < oo, 

T 

||m — Ukh\\L^{I-LP{n)) < Cln — (||u — 7rfeu||is(/.ip(j7)) + \\PhU — m||l'=(/;Lp(17)) + \\RhU — m||l»(/;Lp(i7))) , 

where the projection tt^ is defined in (SB, Ph ■ L^{[2) —^ Vh is the orthogonal L^ projection and Rh '■ Hq([2) —^ 
Vh is the Ritz projection. 


Proof Put e:= u- Ukh = {u - PhT^ku) + (PhT^kU - Ukh) '■= Vkh + ^kh- For 1 < s,p < oo, we have 
\\e\\L‘iI-LP(0)) = sup (e,'i/')7xn, - + ^ = 1, i + ^ = l. 

,peL^'(i;LP'(n)) S s' p p' 

For each such ip, consider a dual problem 

B{pkh,Zkh) = {'Pkh,tp)lxO- 


Thus, we have 


{e,'tp)ixn = {Vkh,'f)ixn + {£.kh,'f)ixn ■= Ji + J2- 
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Using the Holder inequality, the triangle inequality, the stability of the projection Ph in and the approx¬ 

imation properties of and Ph, we find 


Ji < C\\r]kh\\L‘‘ii-,L!‘{n))\\i’\\L‘'{i;Lp'in)) < C\\f]kh\\L‘ii-,LP{n)) = C\\u - Phnku\\L-(i-,LPiO)) 

< c\\u - Phu\\L-{I-,LP{ 0 )) + C\\Ph{u - TTkU^L-il^LPin)) 

< C\\u — Phu\\L‘{i-Lp(n)) + C\\u — TTku\\L‘(i-LP{n))- 

On the other hand, using that B{u — Ukh,Xkh) = 0 for any Xkh € and the properties of the projection 

and the properties of tt/c, we obtain 


M 


J2 — ^{^kh, — ^{Vkh, ^kh) — ^ {Vkh, djZkh) ^ {XTjkh Zkh^ Irn>^ ^ iVkh,m^ \-^kh\m) O 

m—1 

M 

= ^ (m - T^kU,dtZkh)im^n - {X'qkh,Xzkh)in,xQ + (u“ - (tTfeu)", [zkh\m)n 


= -(V(u - Ph7rku),Vzkh)ixO- 

where we used that the first sum vanishes due to (141 al l and the sum involving jumps due to (I41bl l. Using the 
properties of the Ritz projection, integrating by parts in space, and using the Holder inequality and Theorem|6] we 
obtain 


J 2 = -{X{u- PhTrku),Vzkh)ixQ = -{X{RhU- PhTTku),Vzkh)ixn = [RhU - PhT^kU, AhZkh)ixn 

£ C\\Ph{RhU — 

T 

< C\n-\\RhU - T:ku\\L‘-(i-LP{n))H\\L‘'(I-LP'(Q)) 

T 

< C'ln — {\\RhU — u\\L‘{I-,LP{f2)) + ||m — 7rfeu||L = (7;LP(J7))) ■ 

Combining the estimates for Ji and J 2 we obtain the result. 

Corollary4 If the solution u to ^ satisfies u G and fl such that elliptic 

pp 2 ,p _ fggiiidfify holds, then there exists a constant C independent of k and h such that 

T 

||u — UkhWhoil-LP^Cl)) < Cln — (fc'*'''^||M||w'!+i-“(/;LP(r2)) + ^^~^^\W\\L‘(I-WP+^’P(n))) : 1 < s,p < 00 . 

6 Fully discrete results in general norms 

For the future references we provide discrete maximal parabolic regularity results in general norms. For exam¬ 
ple, we use these results to establish pointwise best approximation estimates in lIZTll for fully discrete Galerkin 
solutions. 

Let 17 be a Lipschitz domain and let T = {r} be an arbitrary partition of 17 into cells r (triangles, tetrahedrons, 
quads, or hexahedrons, not necessary quasi-uniform). Let I 4 , be the set of all functions in i7g (17) that belong to a 
certain polynomial space (i.e., Pr or Qr) on each r. As before, we define a fully discrete solution Ukh G X^’J^ by 

B{ukh,I>kh) = {f,‘Pkh)ixn + {uo,p'lh)n for all ipkh G X^’l, (47) 

where 

Xk,I = {vkh ■■ Vkh\i„ GVqiVh), m = 1,2,.. .,M}, for some g > 0, r > 1. (48) 

As in the previous section, we introduce the discrete Laplace operator Ah ■ Vh ^ Vh by 

{-AhVh, x)n = iXvh,Xx)n, Vx G 14, 

and the orthogonal projection Ph - Lf{fl) -G I 4 by 

{PhV,x)f 2 = {v,x)q, VxG14. 
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Let III'III be a norm on Vh such that for some 7 S (0, 7 ) the following resolvent estimate holds, 

|||(z + Z\?j)“^x||| < ^lllxlll. forzGC\r^, (49) 

for all X € 14 , where is dehned in @ and the constant Mh is independent of z. 

For quasi-uniform meshes, this assumption is fulhlled for IH'IH = ||•||LP(^ 7 ) with a constant Mh < C indepen¬ 
dent of h, see as discussed and exploited above. For a weighted norm IH'IH = ||(t'^ •||L 2 (i 7 ) with the weight 
f^xoix) = \/\x — xol"^ + and Mh < (7|ln/i| we established this estimate in lIZTi . and used the corresponding 
result to obtain interior (local) pointwise estimates. Moreover, the resolvent estimate (l49l l is known also to hold in 
LP{n) norms on a class of non quasi-uniform meshes as well, see jj)- 


6.1 Smoothing estimates for the homogeneous problem in general norms 

For the homogeneous heat equation ([T]i, i.e. f — 0 and its discrete approximation Ukh G ^k'h defined by 

B{ukhT ^kh) = i'^0: ^kh,o) '^^kh G (50) 

we have the following smoothing result. 

Theorem 13 (Fully discrete smoothing estimate in general norms) Let |||'||| be a norm on 14 fulfilling the 
resolvent estimate dull. Let Ukh be the solution of dsnii. Then, there exists a constant C independent of k and h 
such that 

sup \\\dtUkhit)\\\ + sup \\\AhUkhit)\\\ + kf,^\\\[ukh]m-i\\\ < ^T^III-P^iUolll, 

tg/m telm ^rn 

for m = 1,2,..., M, where Ph '■ L'^{S1) 14 is the orthogonal L^ projection. For m = 1 the jump term is 

understood as [ufc/i]o = ~ BhUo- 


6.2 Discrete maximal parabolic estimates for the inhomogeneous problem in general norms 

Now, we consider the inhomogeneous heat equation ([T]i, with uq = 0 and its discrete approximation Ukh G 
dehned by 

P{pkh-i^kh) — (,f 1 Tkhf '^Tkh G ^k,h' (51) 


Theorem 14 (Discrete maximal parabolic regularity in general norms) Let IH'IH be a norm on 14 fulfilling 
the resolvent estimate (ED) and let \ < s < oo. Let Ukh be a solution of dlB. Then, there exists a constant C 
independent of k and h such that 



IWAhUkhitWdt 


M 


+ I b^kh] m—1 I 


\m—l 


T 


<CMhlnj[ I WlPhimi dt] , 


where Ph : L‘^[T2) —^ 14 is the orthogonal L^ projection and with obvious notation change in the case of s = oo. 
For m = 1 the jump term is understood as [ukh]o = u^h O' 

The proofs of the above two results are identical to the proofs of the corresponding time discrete results from 
SectionUl provided the resolvent estimate ( l49l l holds. 
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